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We study the electronic states of graphene in piecewise constant potentials using the continuum 
Dirac equation appropriate at low energies, and a transfer matrix method. For superlattice poten- 
tials, we identify patterns of induced Dirac points which are present throughout the band structure, 
and verify for the special case of a particle-hole symmetric potential their presence at zero energy. 
We also consider the cases of a single trench and a p-n junction embedded in neutral graphene, 
which are shown to support confined states. An analysis of conductance across these structures 
demonstrates that these confined states create quantum interference effects which evidence their 
presence. 
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I. INTRODUCTION 

Graphene, a single layer of carbon atoms laid out in 
a honeycomb lattice, is one of the most interesting elec- 
tronic systems to become available the last few years^* 2 -. 
It differs from conventional two dimensional electron gas 
(2DEG) systems in that the low energy physics is gov- 
erned by a Dirac Hamiltonian rather than the more com- 
mon form for semiconductors, characterized by an effec- 
tive mass and a band gap. Because of its unusual prop- 
erties, the possibility of exploiting it for nanoscale appli- 
cations has become an area of intense investigation. 

Nanophysics in graphene may be investigated in sys- 
tems with very small physical size scales, such as quan- 
tum dots 3 or quantum rings 4 9 . More commonly, ma- 
nipulation at the nanoscale is accomplished by the ap- 
plication of electric fields via gate geometries, subject- 
ing the system to potentials varying on a very short 
length scale. Much recent work has focused on p-n 
junctionsiSr-i 4 and p-n-p junctions 15,16 , whose behavior 
in graphene is distinct from that of conventional 2DEG's 
because of the absence of a gap in the spectrum. Very re- 
cently, studies of graphene in superlattice potentials have 
demonstrated the possibility of "band structure engineer- 
ing" of graphen o 17 i 18 . In particular, for one-dimensional 
superlattice a 19 i 2Q , effective potentials may be designed by 
which the number of Dirac points at the Fermi energy can 
be artificially manipulated^i^^"— . 

In this work, we study graphene systems in the pres- 
ence of one-dimensional potentials that are piecewise con- 
stant. Our approach models graphene in terms of the 
Dirac equation, which is valid for potentials that vary 
slowly on the lattice scale, and for which inter valley scat- 
tering may be ignored. The allowed energy states for such 
systems may be found by a transfer matrix method de- 
scribed below. For a one-dimensional superlattice, such 
problems are reminiscent of dynamical systems with pe- 
riodically varying parameters. We find that the resulting 
band structure generically supports a collection of Dirac 



points at a variety of energies. Dirac points generally ap- 
pear at energies satisfying e = Ed/hv F = nir (E = elec- 
tron energy, d = superlattice period, v F = Fermi veloc- 
ity) for n an integer, at q = 0, where q is the momentum 
perpendicular to the superlattice axis. For particle-hole 
symmetric potentials, several Dirac points may appear 
at e = for non- vanishing values of g2i ~ 23 i 26 ~ ». Generi- 
cally, other Dirac points are present for e, q ^ 0. 

We also study geometries with a single p or n region, 
or with a p-n junction region, embedded in a neutral 
graphene background. We refer to these as a single trench 
geometry and an embedded p-n junction, respectively. 
The trench geometry supports localized solutions, at en- 
ergies for which the electron states are evanescent in the 
neutral graphene for a given k y ^ 0. These states dis- 
perse as a function of k y such that a finite number of 
them crosses zero energy. The embedded p-n junction 
supports localized solutions as well, which may be inter- 
preted in terms of a pair of overlaid trench spectra. If the 
two trenches are particle- hole conjugates, the individual 
trench spectra cross zero energy at the same value of k y ; 
however, because of tunneling between the trenches these 
become avoided crossings rather than Dirac points. At 
large k y the gaps at the avoided crossings may be quite 
small. 

Finally we consider transport across these systems 
within the Landauer formalism, for which the transfer 
matrix method is easily adapted to yield the conduc- 
tance of the system. We demonstrate that the conduc- 
tance has distinctive signatures indicative of the confined 
states which may be contained at non- vanishing k y for a 
given embedded structure. 

This article is organized as follows. In Section [II] we 
express the problem of a Dirac particle in a periodic po- 
tential in terms of a path-ordered product, and find the 
band structure of the simplest case of alternating positive 
and negative potential regions. Section IIIII develops the 
transfer matrix method for handling more general piece- 
wise constant potentials. In Section IIVI we apply the 
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method to compute the spectrum of localized solutions 
in trench and embedded p-n junction geometries. Section 
W\ applies the results of the transfer matrix calculations 
to compute the conductance of these system. Finally, we 
conclude with a summary in Section IVII 



II. PERIODIC POTENTIAL 

In the vicinity of a Dirac point at wave vector K, we 
write the wave vector k = K + Afc, and the two compo- 
nent wavefunction as ij) = e lK x (™). Assuming a scalar 
potential V(x) which is independent of the coordinate y, 
the wavefunction <j) = (™) may be chosen an eigenstate 

of the Hamiltonian H = e - iK ' x H e iK ' x : 



variety of problems involving linear equations with peri- 
odic coefficients 24 , such as the classical mechanics setting 
of a pendulum with periodically varying length. 

For a continuous periodic function V(x) the path or- 
dered exponential is not analytically tractable. However, 
many of the essential features are reproduced by consid- 
ering a piecewise constant V(x) of the form 



V(x) 



Vn 



if 
if 



< x < d 1 
d>i < x < d x 



d 



(8) 
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We then have A = A 2 A l5 where A ■ = exp(M-d^), where 
the matrices Mj (j = 1,2) are as in Eq. [2j with Kj = 
(E — Vj)/hv F in the off-diagonal elements. The period of 
V(x) is d 1 + d 2 . We then have 



H 



V(x) hv F (-id x + iq) 

hv F (-id x -iq) V{x) 



(i) 



where q = Afc . The eigenvalue equation for <\> may then 
be recast in matrix form as d x (j) = M{x) (/>, where 



M(x) 



and 



k(x) 



-q zk(x) 
in(x) q 



E - V{x) 



(2) 



(3) 



The solution is 

^u(x)\ __ 
v(x) 



Pexp dx' M{x') 




where P denotes the path ordering operator, which places 
smaller values of x to the right. 

We presume V(x) = V{x + d) is periodic, in which case 
the Bloch condition becomes 



u(d) 
v(d) 



u(0) 
v(0) 



(5) 



where 6 X is the Bloch phase. This leads to the condition 
det[A-e^] = , (6) 

where 



/ d 

{h 

\ n 



A = Pexp dxM(x) 



(7) 



Note that Aft -a x Ma x , where a x is the Pauli matrix, 
and that TrM = 0. Since det(A) = 1, the eigenvalues 
of A are given by ( = \T ± VT 2 - 4, where T = Tr (A). 
The energy bands, where the Bloch condition can be sat- 
isfied, correspond to the regime \T\ < 2. For \T\ > 2 the 
wavefunctions grow exponentially and are not normaliz- 
able. The mathematical structure here is common to a 



sin a a 

Aj = cos otj'I-\ • Mjdj 



(9) 



with a - — dj 



k?- — q 2 . The resulting trace is 



smai sinfto 



T = 2 cos a-, cos a 9 + 2 d, d 9 (q 2 — fti k 9 ) 

1 1 1 1 v 1 l) a r a 2 

(10) 

The Bloch condition is satisfied when \T\ < 2. It is con- 
venient to define the dimensionless quantities 



-> = Uv, - v 2 ) ■ 



d 1 + d 2 



(ii) 

(12) 



( 4 ) where (V) = (d x Vi + d 2 V 2 )/(d 1 + d 2 ) is the average 



potential. We also define the dimensionless wavevec- 
tor = q(d 1 + d 2 ). When V is constant, one finds 



T = 2 cos 
fled provided \E 



6 2 

u y ■ 



and the Bloch condition is satis- 



V\ > hv F \q\. This, of course, is the 
condition that the energy E lies within the Dirac cone. 

Next, consider V x ^ V 2 , but with q = 0. The trace is 
then T = 2 cose, and there are allowed states at all en- 
ergies. The Bloch condition is satisfied marginally when 
e = n7r, corresponding to the touching of two energy 
bands which disperse with g, i.e. to the vanishing of a 
band gap as q — ^ 0. Indeed this result holds for arbi- 
trary V(x), since M(x) = n{x) a x in this case, hence 
[M(x), M(x')\ = and the path ordering is superfluous. 

Now consider the case d x = d 2 . The spectrum then 
exhibits a particle- hole symmetry relative to E = (V). 
When E = (V), i.e. e = 0, we have 



T = 2 



26 y sin \ 



2 



01 



Clearly T > 2, with the marginal condition T 
taining when 



<9 = ±^uj 2 - 47r 2 n 2 



(13) 



2 per- 



(14) 
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FIG. 1: (color online) Contour plot of T = Tr (A) for d\ = d 2 
and uj = ^tt. Red regions correspond to T > 2, and olive 
regions to T < — 2. The energy bands are identified with the 
blue regions, where \T\ < 2. 



where n is an integer, as well as a solution at = 0. 
Since q is real, there is a maximum allowed value for n: 



FIG. 2: (color online) Contour plot of T ■ 
and d\ — | d 2 . 



■ Tr (A) for uj = fir 



where a = (e — ^ ) 2 — Oy . This limit is distinguished 

by the fact that all the Dirac points of T(e,0 y ) lie along 
the line = 0, as shown in Fig. 2J 



UJ 

2tt 



(15) 



where the brackets indicate the greatest integer function. 
Thus, two new Dirac points open up each time uj in- 
creases by 27r, and there are 2n max + l Dirac points, in- 
cluding the one at y = 0. Similar results have previously 
been obtained by other method a 21 i 22 . 

In Figs. Q] and [2] we plot the contours of the function 
T(e,O y ) for uj = ^-tt for d 1 — d 2 and d 1 = \d< 1 , respec- 
tively. Note the particle hole symmetry present in the 
former case, which is broken when d x ■=£ d 2 . Note also 
the band touchings at e = nir for = in both cases. 
In Fig. [3] we again plot the band structure as in Fig. [TJ 
but showing the high energy structure as well. At en- 
ergy scales \e\ ^> 1, the spectrum may be understood as 
the intersection of the two Dirac cones \e 
\s-uj\> y . 

As a limiting case of the model, consider d 2 
V 2 d 2 = —hv F Q finite. The potential is then 



uj\ > y and 



with 



III. TRANSFER MATRICES AND LANDAUER 
CONDUCTANCE THROUGH POTENTIAL 
SLABS 

Consider a piecewise constant potential V(x) as de- 
picted in Fig. [5] in which V = Vj over a slab of width d^ 
with j = 1, . . . , N. We presume the slabs are bounded 
by pristine regions of graphene with V = 0, and we label 
these outer regions and N 
of the slabs may be defined by the relations 



1. Wavevectors kj for each 



E 

hv 



(18) 



where each k- > is nonnegative. Here, a- = ±1 is a 
binary variable in each slab which indicates whether the 
energy E is above (n-type) or below (p-type) the Dirac 
point 



which occurs at E = We take a = cr N+1 and 



AT+l 



0. 



In each region j, the solution to the Dirac equation, 



V (x) = n -^h- J2 6 {h n )\> (16) 



n= — oo 



i.e. a Dirac comb with a uniform compensating back- 
ground. The trace in this case is given by 



T = 2 cos Q cos a — 2 sin Q sin a • 



£ — fir 



a 



M, (17) 



iHv F a • + VA) = Eip 



(19) 



is either propagating or evanescent. Due to transla- 
tion invariance in the ^/-direction, we are free to fix the 
wavevector q along this axis; all solutions include a fac- 
tor e iqy . The local Fermi wavevector fc - in each region is 
related to the local density by rtj = /c 2 /47r, and the sign 
of the carrier charge is —a-. 
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FIG. 3: (color online) Contour plot of T = Tr (A) for u = ^fir FIG. 4: (color online) Contour plot of T = TV (A) in the 



and di — cfe 



d2 —> limit with = —V2d2/hv F = Itt held constant 



If fc ■ > g, the wavefunction in slab j propagates as a exponential) behavior along the x-axis: 



plane wave along the x-direction as well. The solution is 
of the form 



V2 



e ik jX cos e J + 



(20) 



e -ik jX cos 9 j I e ik jV sin Q. 



where 6 G [ — § , f]. Matching the ^/-dependence gives 
us the condition 



q = k j sin ■ if k, > q . 



(21) 



To solve for the wavefunction everywhere we must 
match the value of each component of ip(x) at the bound- 
aries between slabs. Accordingly, from Eq. [20j we define 
the matrix 



Ma 



1 



TS^ 6 3 



(22) 



cos(^ 
iaj sin ( 



v 3 - 



e k j xcot Vj + 



(24) 



(?) 



■kjX cot ifj I ik-ycsc^- 



where cpj G [ — § 5 §]• Momentum conservation along y 
requires 



# = kj esc if kj < q . 
For the evanescent slabs, we define 

f cos(^) sin(^) > 
\icj j sin (^) io.caatff)) 



(25) 



(26) 



and free propagation across an empty slab is described 

by 



'e k j d j cot 











(27) 



which we shall invoke presently. Free propagation across 
a slab is described by the transfer matrix 



, gikjdj cos 6 



o 



— ikjdj cos 6 



(23) 



If k ■ < q, the solution in slab j has an evanescent (pure 



The equations guaranteeing continuity of the wave- 
function at the slab boundaries are 



M 



j-i 



M, 



(28) 



Here, the superscripts R and L refer to the right and left 
boundaries of the slab. They are related, within a given 
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FIG. 5: Transport through a series of potential slabs. The 
potential V(x) is piecewise constant, given by Vj in the j th 
slab. 



slab, by the free propagation transfer matrix, 




(29) 



We can now compute the transfer matrix M for a gen- 
eral configuration of slabs, which is defined by the rela- 
tion 



f A L 



M 



N-\-l / 



^0 



OR 

^0 ' 



(30) 



In terms of the previously defined quantities, the total 
transfer matrix is 

M = -Wjv+i {M N K N M^ 1 ) ■■■{M 1 K x M^)M . (31) 

The recurring matrices Qj = MjKjM^ 1 may be parame- 
terized by 6j and Xj = kdcosOj for a propagating region, 
and by (fj and Uj — kd cot pj for an evanescent region, as 
detailed in the Appendix. As is well-known, the transfer 
matrix relates data on the left to data on the right of a 
given region. Incoming and outgoing flux amplitudes are 
related by the ^-matrix, 



(32) 




Appealing to the unitarity of S, we then have 

M n =t*~ 1 
M 21 



M 12 =r't'~ 1 



-t'- 1 r 



M 



22 



t' 



(33) 
(34) 



As discussed in more detail below, the dimensionless 2- 
terminal Landauer conductance, in units of e 2 /h per spin 
channel, is given by 



> 
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FIG. 6: Energy bands for a single trench of width d and 
potential V = 14 hv F /d. The hatched region depicts the Dirac 
cones in the bulk. Solid lines show dispersions for bands of 
localized states in the gap between the upper and lower Dirac 
cones. 



IV. LOCALIZED SOLUTIONS OF EMBEDDED 
GEOMETRIES 

We first consider situations where the outer regions 
contain evanescent waves, so the wavefunction is lo- 
calized within the slabs. In this case, we must have 
B = and A N+1 = 0, if cotcp > 0, which requires 
Mn = 0; for the case cot p < 0, we must have A = 
and B N+ i = 0, leading to M 22 = 0- This condi- 
tion defines the one-dimensional dispersion E(q), where 
q = k esc ip 



kj sin t 



A. Single trench 

For the case where there is a single trench of width <i, 
the condition Mn = is 

cos ip cos 1 cos Xi + ( sm #i — G \ sm Po) sm Xi = 

(36) 

or 

cos(/? cos(/? 1 coshcj 1 + (l— a o x sin^ sin^) sinhcj 1 = . 

(37) 

The condition M 2 2 = is obtained by reversing the sign 
of Xi or ^i- We now define the dimensionless poten- 
tial v = Vd/hv F , energy e = a k d = v + a^d, and 
wavevector = qd. The eigenvalue condition is then 
F(e, v, 0y) = 0, where 



1*1- 



(35) 



2 - e 2 cos Xi 
(38) 
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FIG. 7: Energy bands for a p-n junction embedded in neutral 
graphene. Both p and n regions have width d, with Vi = 
—Vi — 14:hv F /d. The hatched region depicts the Dirac cones 
in the bulk. Solid lines show dispersions for bands of localized 
states in the gap between the upper and lower Dirac cones. 



and 



Xi 



e 2 



(39) 



This form for F holds for either sign of q and is extended 
to the regime (e — v) 2 < 2 via analytic continuation. If 
we search for eigenvalues at e = 0, we obtain the condi- 
tion 



Xi 



(y 2 - x?) 1/2 



■tanxi , 



(40) 



a graphical analysis of which shows a new root emerging 
each time ud passes through nir. Thus, there are [z//7r] 
solutions with e = which lie at finite values of . In 
Fig. [6j for example, we have v = 14 and thus there are 
[14/7r] = 4 bands which cross e = at finite . 



B. Embedded junction 

Consider next the case of an embedded junction with 
a double step, 



V[x) 



if x < 

Vi if < x < d t 

V 2 if d x < x < d t + d 2 

if d x + d 2 < x . 



(41) 



When V x and V 2 are of opposite sign, this corresponds 
to a p-n junction embedded in neutral graphene. We 



FIG. 8: Low energy detail for the embedded p-n junction from 
Fig. [71 showing that E = always lies within a gap between 
subbands. 



define e = cr k (d 1 + d 2 ), y = q(d t + d 2 ), and v- = 
Vj(d 1 +d 2 )/ hv F . We furthermore define x- = dj/(d 1 J rd 2 ), 
so x 1 + x 2 = 1. Applying the transfer matrix formal- 
ism, we once again arrive at an eigenvalue condition 
F(e, v 1 ,v 2 ,x 1 , y ) = 0, this time with 



F(e, v 1 ,v 2 ,x 1 ,0 y ) = x 1 (el -eie-vS) 



smxi 
Xi 



J x 2 



fn2 ( \\ sin X2 

+ X 2 [°y ~ e ( £ ~ V 2)) * COS Xl 

^1^2 (fy ~ (e-^l)(e-^2)) 



x 2 

sinxi sinx 2 



cosxi cos X 2 



where 



Xi 



6> 2 



(42) 



(43) 



and where we continue to express energies and wavevec- 
tors as dimensionless quantities. As with the single 
trench case, the expression for F is analytically continued 
to the regimes where (e — z^) 2 < 2 . 

Generally speaking, when two semi-infinite undoped 
regions are separated by a sequence of doped slabs, in 
addition to the Dirac cones in the bulk one obtains local- 
ized subbands with energies lying within the gap between 
the upper and lower Dirac cones. The wavefunctions for 
these subband states are localized within the slabs. In 
Fig. [6j we plot the Dirac cones and subband dispersions 
for the case of a single strip of width d with potential 
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FIG. 9: The bulk Dirac cones for V — are hatched in black. 
Localized states in the p trench (Vi = 12hv F /d) are present 
in regions I and III. Localized states in the n trench (V2 = 
—7hv F /d) are present in regions I and II. 



FIG. 10: Subband structure and bulk Dirac cones for the 
p-n junction with d\ = cfe = d and V\ — 12hv F /d, V2 = 
—7hv F /d, shown in black. Superimposed in blue and red are 
the subband structures for the individual Vi and V2 barriers, 
respectively. 



V = 14hv F /d. In Figs. [7]and[8j we investigate the case 
of a p-n junction embedded in neutral graphene, where 
both p and n regions have the same width d, and again 
we choose V = ±14hv F /d for purposes of illustration. 
Apparently E = always lies within a gap between sub- 
bands. Note that the parameters of the p-n junction are 
such that the spectum is particle- hole symmetric. 

The subband structure of a general p-n junction can 
be understood in terms of the simple diagram in Fig. 
[HI In the figure, we sketch the three Dirac cones for 
the undoped bulk region (black lines), the p-type slab 
(blue lines, V 1 d 1 /hv F = 12), and the n-type slab (red 
lines, V 2 d 2 /hv F = — 7). For the purpose of generality, we 
consider a system which has no particle-hole symmetry. 
Outside the black hatching, there are no bulk states. In 
regions I and III, there are localized states living along 
the p-trench. Similarly, in regions I and II, there are lo- 
calized states living along the n-trench. In region I, there 
can be resonant tunneling between the p and n trenches, 
and the localized states live in both trenches. In region 
IV, which lies outside all three Dirac cones, there can 
be no states. In Fig. [lOl we plot the subband struc- 
ture for this p-n junction, with equal thicknesses of the p 
and n regions, superimposing on those curves the results 
for the individual single barrier problems. A comparison 
with the sketch in Fig. [9] reveals the basic physics of the 
energy spectrum. One may see the confined states associ- 
ated with each of the two trenches, one set dispersing up 
and the other down. These confined states admix where 
they are degenerate for the individual trenches, leading 
to anticrossings. In the case where the two trenches are 
equal and opposite (e.g., Figs. [3 and [8]), this leads to a 



set of anticrossings around zero energy, rather than the 
extra Dirac points one finds in the superlattice case. 



V. TWO-PROBE LANDAUER CONDUCTANCE 

In this section, we consider transmission through a slab 
geometry for the more general case where the graphene 
leads are doped. Consider the geometry of Fig. [5] De- 
fine a vertical (y-directed) surface E somewhere to the 
right of the N th slab. The current across that surface 
will be a sum of three contributions. First, there will be 
carriers injected from the left side, with chemical poten- 
tial /i L , moving at angle 6 G [ — ^tt, \tt\ , which make 
it past the slabs with probability T(E,0). These are 
all right-movers. Second, there will be carriers from the 
right lead, with chemical potential /i R , moving at an an- 
gle G [^7r, |7r] . Because of this range of all these are 
left-movers. Finally, with probability R'(E,0), each of 
these left-moving carriers reflects off the barrier region, 
scattering into a state with 0' = ir — 0. Putting this all 
together, 



00 tt/2 

J E = L E n E • j, dE Af(E) J ^ {t(£, 0) J(E, 0) f(E - h ) 

-00 -7T/2 

J(E, 7T + 6) + R(E, 7T + 6) J(E, -6)] f(E- h )j, 

(44) 
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FIG. 11: Average transmission coefficient across a single 
trench of width d\ — 100 nm as a function of density m (in 
units of 10 12 cm -2 . From top to bottom: no = 1, 2, 3, 4, 
5, and 20 (xl0 12 cm -2 ). Negative densities correspond to a 
p-type (hole-doped) trench. 



FIG. 12: Average transmission coefficient across an embedded 
p-n junction with d\ — d<z — 50 nm as a function of density 
n\ — —ri2 (in units of 10 12 cm -2 . From top to bottom: no — 
1, 2, 3, 4, 5, and 20 (xl0 12 cm- 2 ). 



where L s is the length of the surface and Af(E) is the 
density of states, 

WE > - SiKJ ' < 45 » 

and g = 4 accounts for spin and valley degeneracy. 
J(E, 9) is given by 

J(E, 9) = -ev F sgn(E) (x cos 9 + y sin 0) . (46) 

We take n E = x. Unitarity of the ^-matrix yields 
T(E, 9) + R(E, 7T + 9) = 1. Expanding in the difference 
Ml — Mr, we obtain at T = 

I s = ^V{T(E r) )^L,, (47) 

where V = (/i L — /i R )/e is the voltage drop. The average 
transmission coefficient is defined as 

tt/2 

(T(E))= Jd9 T(E,9) cos<9 . (48) 
o 

The carrier density at zero temperature is n = gkl/Air, 
so the Fermi energy is E F = hv F k F = hv F ^ Ann/ g. Thus, 

I(V) = ^ (T(E F )) (47T5U) 1 / 2 L y V . (49) 
Defining the effective channel number, 

N c = ^k F L y = ^gn/^fl 2 L y , (50) 



the two-terminal Landauer conductance becomes 

G=^=N c x^-(T(E F )). (51) 
Expressed in a convenient set of units, the resistance is 

= 0.114kfi • (T(£ F )) -1 • n~ 1/2 ■ I' 1 , (52) 

where n is the dimensionless carrier density (in the leads) 
in units of 10 12 cm -2 and L y is the length of the barrier 
region in units of microns. 

For a single trench of width d x and carrier density n 1 
embedded between graphene leads of carrier density n , 
we have M = M _1 Q 1 M , where Qj Mj K-MJ 1 , and 

Re.M n = cosher (53) 
I m Mn = a o cr i sinh^ tan cp 1 sec 9 (54) 
— sinh uj 1 sec cp 1 tan 9 . 

Here we have assumed fc 1 < q so the trench contains 
evanescent solutions, but by a simple analytic continua- 
tion our result also holds for the k x > q case. 

For an embedded junction consisting of two consecu- 
tive slabs (n 1 ,(i 1 ) and (n 2 ,d 2 ), the total transfer matrix 
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is M = Mq 1 Q x Q 2 M , resulting in 
Re jM n = cosh u 1 cosh uo 2 (55) 
+ sinh lj 1 sinh uo 2 ( sec (p ± sec tp 2 — cr 1 a 2 tan tp x tan y? 2 ) 
ImA^! = <r cr 1 sec # sinnu^ tan y? 1 cosh uj 2 (56) 

+ <j <r 2 sec #o cosh uj 1 sinh uj 2 tan <^ 2 
— tan #o ( cosh uj 1 sinh uj 2 sec y? 2 + sec (p 1 sec ^ cosh co> 2 ) . 

In Fig. [11] we plot the average transmission coeffi- 
cient (T(E F )) as a function of the density n 1 for a sin- 
gle trench with d x = 100 nm, for several different values 
of the lead density n . There are three features wor- 
thy of note. First, there is a nonzero transmission even 
when n 1 = 0, due to the presence of evanescent solu- 
tions within the trench 25 . Second, there is (obviously) 
perfect transmission when n 1 = n . Third, the trans- 
mission falls off rapidly as n 1 — >• 0, and is smooth for 
< n 1 < n , which is a regime where the carriers in 
both regions are of the same sign, yet total internal re- 
flection can occur at the interface, with the critical angle 

given by sin#Q = r^/no . Finally, when n 1 > n 
or when the carriers in the leads and the trench are of 
opposite sign, one observes quantum oscillations in the 
transmission, due to constructive and destructive inter- 
ference effects from scattering at the walls of the trench. 

It is interesting to compare Fig. [11] with results in Ref. 
HTL which reports the measured conductance through 
an electrostatically defined trench potential in a wide 
graphene ribbon as a function of the gate voltage defining 
the trench. As in our results, the experiments indicate a 
rapidly falling conductance when n 1 has the same sign as 
n , as n 1 approaches zero, and after passing through zero 
one observes oscillations due to quantum interference. A 
prominent difference is the absence in the experimental 
results of a rising background around which these oscilla- 
tions occur. In our model, this appears because the den- 
sity of states in the trench increases as \n 1 1 increases. It is 
likely that the difference between our theoretical results 
and those observed in Ref. [ll| are due to the relatively 
smooth potential induced by the remote gates in the ex- 
periment. In particular this can yield relatively large 
"depletion" regions between the trench and the graphene 
leads which may dominate the resistance of the device, 
and overwhelm the increase in conductance related to the 
trench density of states. In principle, one could model the 
device as a series of piecewise constant potentials and use 
the transfer matrix method we have developed to obtain 
a more quantitatively accurate model. 

In Fig. [12] we plot (T(E F )} as a function of n 1 for 
an embedded p-n junction with d x = d 2 = 50 nm and 
n 2 = — n 1 . For such a symmetric junction the average 
transmission is symmetric under n 1 — >• —n 1 . The step- 
like features present in the transmission here represent 
scattering effects due to the localized solutions and the 
avoided crossings in their dispersions, as discussed in the 
previous section. Thus the bound states support an (in 
principle) measurable signature in conductance through 
the system. 



VI. SUMMARY 

In this paper we have studied the spectra of and trans- 
port properties through graphene in the presence of sev- 
eral piecewise constant potentials. We demonstrated 
that these may be investigated in a straightforward man- 
ner using a transfer matrix approach. For the case of 
super lattice potentials, we demonstrated that exact so- 
lutions of the Dirac equation generically contain Dirac 
points centered at zero energy for particle-hole symmetric 
potential ^ 21 ! 22 , and that these are special cases of Dirac 
points that appear throughout in the band structure of 
superlattice potentials. We also considered the cases of 
a single trench and an embedded p-n junction in neu- 
tral graphene, both of which may support confined states 
at large enough values of k y , the momentum along the 
translationally invariant direction. An analysis of con- 
ductance through these structures shows quantum inter- 
ference effects which may be interpreted as signatures 
of these bound states, in principle allowing these unique 
features of graphene nanostructures to be seen experi- 
mentally. 

Note added: While this manuscript was in the fi- 
nal stages of preparation, a manuscript by Barbier, 
Vasilopoulos, and Peeters appeared 2 ^ which also consid- 
ers the effect of a superlattice on monolayer graphene 
and the emergence of additional Dirac points. We are 
grateful to F. M. Peeters for subsequently alerting us to 
the substantial previous work which he and his collabo- 
rators have contributed on the subject of Dirac subband 
spectra and transport in the presence of one-dimensional 
potentials^"—. We are also grateful to Y. P. Bliokh for 
drawing our attention to Ref. 29, which also addresses 
subband formation in the presence of a superlattice struc- 
ture. 



VII. APPENDIX 

We define Q = MKM~ X for a given slab. Within a 
propagating slab, we have 

n _ fcosx + sin x tan # ias'mxsecO \ , . 
\ icr sin x sec cos x — sin x tan 6 J ' ^ ' 

with x — kd cos 9, while for a nonpropagating slab, we 
have 

n _ /cosh uj + sinh uj sec ip ia sinh uj tan ip \ 
y ia sinh uj tan cp cosh u — sinh uj sec cpj 

(58) 

with uj = kd cot ip. These expressions are useful when 
constructing the transfer matrix for a series of slabs, 
which is given in Eq. [31] 
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